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1. Introduction 

In this paper we present a new determination of the QSO Luminosity Func- 
tion (LF) and its evolution, based on new data coming from the Homoge- 
neous Bright QSO Survey (HBQS, partly published in Cristiani et al. 1995). 
We also take into account the most relevant already existing surveys. 

The HBQS was started in 1989 in the framework of an ESO Key- 
programme. The survey covers a total of 555 deg 2 subdivided in 22 ESO or 
UKST fields at high galactic latitude around the south galactic pole. Two 
Schmidt plates for each bandpass U, B' or Bj, V , R or OR and / have 
been obtained at the UKST or ESO Schmidt telescope, usually within a few 
months interval, in order to minimize the effects of variability. The plate 
material has been scanned on the COSMOS microdensitometer. According 
to the ESO/UKSTU numeration, the 22 fields which have been observed 
are: 287, 290, 295, 296, 297, 349, 351, 355, 406, 407, 408, 410, 411, 413, 468, 
469, 470, 474, 479, 534, SA94, SGP. 

We have selected as candidates all the UVx "not-extremely-extended" 
objects satisfying a type of modified Braccesi less-restricted two-color cri- 
terion (La Franca, Cristiani & Barbieri 1992). The magnitude intervals in 
which the selection is virtually complete vary from field to field between 
15.0 < B < 17.3 and 15.0 < B < 18.8. The sample extracted from these 
data includes 327 QSOs (M B < -23, q = 0.5, and H = 50 Km/s/Mpc). 
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Figure 1. The QSO number counts. Continuous line: LDLE model C. Dashed line: PLE 
model B. Dotted line: PLE model A (see Table 1). 



The other samples we have used for the computation of the QSO LF 
are: A) The Edinburgh QSO Survey (EQS, Goldschmidt et al. 1992), which 
covers 333 deg 2 in the interval 15.0 < B < 16.5; B) The Large Bright QSO 
Survey (LBQS, Hewett, Foltz and Chaffee 1995), which covers 454 deg 2 
in the interval 16.5 < B < 18.5; C) The SA94 Survey (La Franca et al. 
1992), which covers 10 deg 2 in the interval 15.0 < B < 19.9; D) The 
Durham/AAT Survey (Boyle et al. 1990), which covers 11.9 deg 2 in the 
interval 16.0 < b < 20.9; E) The ESO/AAT (Boyle, Jones and Shanks 
1991), which covers 0.85 deg 2 from J = 18.0 down to J = 21.75 - 22.0; F) 
The (ZM) 2 B Survey (Zitelli et al. 1992), which covers an area of 0.69 deg 2 
down to J = 20.85 and 0.35 deg 2 down to J = 22.0. 

2. The QSO Counts 

The B number-magnitude counts, in the interval 0.3 < z < 2.2, for all the 
samples used in the computation of the optical LF are shown in Fig. 1. 
As first indicated by Goldschmidt et al. (1992), for magnitudes brighter 
than B = 16.4 the QSO surface density turns out to be a factor 2.5 higher 
than what measured by the PG survey (Schmidt & Green 1983). The data 
from the HBQS collect 7 QSOs at magnitudes brighter than B = 16.4, 
corresponding to a surface density of 0.013l;oo5 deg~ 2 . With the addition 
of the EQS the surface density at B < 16.4 becomes 0.016 ± 0.005 deg~ 2 , 
corresponding to a total of 14 QSOs over an area of 888 deg 2 . In a larger 
area (~ 6800 deg 2 ), the PG survey collects 41 QSOs corresponding to a 
QSO surface density of 0.006 ± 0.001 deg' 1 . 
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3. The Luminosity Function 

In the interval 0.3 < z < 2.2 the QSO LF has been usually parameterized 
with a double power law Pure Luminosity Evolution (PLE) model 

$(Mb,z) = 10 0.4[M B -M s ( z )|( a +l)^_ 10 0.4[M s -M s ( z )|(/3+l) ) 

where a and (3 correspond to the faint-end and bright-end slopes of the 
optical LF, respectively. With this parameterization the evolution of the LF 
is uniquely specified by the redshift dependence of the break magnitude, 
Mb{z) = M B — 2.5/clog(l + z) corresponding to a power law evolution in 
the optical luminosity, L* oc (1 + z) k (see Boyle, Shanks, Peterson 1988). 

In order to compute the QSO LF, we have simulated, via Monte Carlo 
techniques, the effects of the various stochastic processes that in reality 
contribute to determine the observed magnitudes and colors. 

The apparent B magnitudes of the QSOs have been computed following 
the K-correction according to Cristiani & Vio (1990), and computing the 
average galactic absorption A B according to the HI maps: 

B = M B + 5 log[A(z)c/H ] + 25 + K(z) + A B + e(B) 
where e(B) simulates a noise with Gaussian distribution with variance a B = 
<7p hot (B) + 0£ i7 - cr B takes into account both the photometric errors of the 
survey a p h t(B) and the apparent magnitude dispersion due to the QSO 
spectral slope dispersion o Brf = 2.5cr 7 log(l + z). 

The apparent colors have been generated following the average QSO 
colors-redshift dependence F\^{z) (Cristiani & Vio 1990). Between each 
couple of bandpasses centered at Ai and A2 the resulting apparent color 
is: Ci 2 = ^1,2(2) + £ Fi2 where £f 12 simulates a noise with Gaussian dis- 
tribution with variance a\ x 2 = a\ + a\ + op a - <J Fi 2 takes into account 
the photometric errors of the two bandpass o\ and 02 respectively, and 
the QSO color dispersion due to the spread a 7 in the slope of the spectra 
<7/? 7 = 2.5<t 7 log The flux limits and selection criteria of the survey have 
been applied to the resulting magnitudes and colors in order to select the 
"observed" QSOs. 

4. Results and Discussion 

A total sample of 1022 QSOs has been used. The fits have been carried out 
by minimizing the x 2 statistics derived from the comparison of the observed 
(B, z) distribution with 2000 simulations of each theoretical LF model ( see 
Table 1). The significance of the fitting has also been tested for goodness- 
of-fit in the interval 0.3 < z < 0.6, in which the fitting probabilities have 
been computed using the 2D-KS test. We first focus our discussion on a 
q = 0.5 universe, we will later see how our results change for q Q = 0.1. 
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Figure 2. The QSO LF. Dotted lines: best fit PLE, model B (see Table 1). 



In Fig. 2 our best fit PLE model (B) is shown. Our PLE model B gives 
a satisfactory fit of the (B, z) distribution with a global x 2 probability 
of 0.21. The PLE model A by Boyle (1992), compared with the present 
data, provides a lower but still acceptable probability of 0.12, however with 
one additional parameter, z cut = 1.9, representing the redshift at which the 
luminosity evolution "switches off" . At redshift greater than z cut no further 
luminosity evolution takes place. Our B model has a flatter bright slope (3 
(—3.7 compared to —3.9) than the A model. The evolution parameter k is 
also smaller (3.26 compared to 3.45). 

This difference in the parameter (3 is originated by a larger density of 
luminous QSOs (expecially at z < 0.6) in comparison with the previous 
data derived from the PG sample (compare Fig. 2 with Fig. 3b of Boyle 
1992). The higher evolutionary rate k determined by Boyle (1992) is a result 
of the introduction of the redshift cut off at z cut = 1.9, which with our data 
results unnecessary. 

Both A and B models provide an inadequate simulation of the distri- 
butions of the observed data for the low redshift domain, 0.3 < z < 0.6. 
For this subsample the 2D KS statistics test rejects model A at the 0.01 
level and model B at 0.02 level. At magnitudes brighter than Mb = —25, 
in the interval 0.3 < z < 0.6, 32 QSOs are observed, while 16 and 19 QSOs 
are predicted by model A (a 4 a discrepancy) and B (a 3 a discrepancy) 
respectively. But as the low redshift (z < 0.6) subsample contain only 5 
per cent of the complete data set, in a global comparison of the whole data 
sample in the interval 0.3 < z < 2.2, the models follow the evolution of the 
larger fraction of QSOs at higher redshift, allowing the x 2 probability of the 
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Model 




fci 


k 2 


a 


P 


M* 


x 2 


2DKS 




mag^Mpc^ 3 












z < 2.2 


z < 0.6 


A 


6.5 x 1(T 7 


3.45 




-1.5 


-3.9 


-22.4 


0.12 


0.01 


B 


1.1 x 1(T 6 


3.26 




-1.39 


-3.72 


-22.3 


0.21 


0.02 


C 


9.8 x 1(T 7 


3.33 


1.2 


-1.45 


-3.76 


-22.3 


0.53 


0.09 


D 


1.0 x 1(T 6 


3.29 




-1.39 


-3.69 


-22.4 


0.14 


0.03 


E 


1.1 x 1(T 6 


3.29 




-1.39 


-3.72 


-22.3 


0.18 


0.02 


F 


1.1 x 1(T 6 


3.27 




-1.38 


-3.72 


-22.3 


0.18 


0.01 


G 


8.5 x 1(T 7 


3.19 




-1.53 


-3.87 


-22.4 


0.13 


0.01 


H 


1.2 x 1(T 6 


3.00 




-1.37 


-3.81 


-22.4 


0.05 


0.01 


I 


4.5 x 1(T 7 


3.55 




-1.6 


-3.8 


-22.6 


0.004 


0.01 


L 


3.4 x 1(T 7 


3.44 




-1.57 


-3.75 


-22.9 


0.36 


0.08 


M 


4.0 x 10~ 7 


3.45 


0.9 


-1.50 


-3.74 


-22.8 


0.49 


0.11 




la errors: 


±0.07 


±0.3 


±0.07 


±0.13 


±0.2 







NOTES: For q = 0.5: A) Boyle 1992; B) PLE; C) LDLE; D) PLE without HBQS; E) 
PLE without EQS; F) PLE without LBQS; G) PLE with cr 7 = 0.3 H) PLE with 
cr 7 = 0.5. For go = 0.1: I) Boyle 1992; L) PLE; M) LDLE. 



global fit to reach satisfactory levels. No significant difference is obtained 
by fitting the data excluding each of the three bright samples (HBQS, EQS 
and LBQS) in turn (models D, E and F). 

It is interesting to see how the inclusion of a spread in the theoreti- 
cal average QSO spectral slope modifies the best fit luminosity function. 
As expected (Giallongo and Vagnetti 1992) larger values of the spread of 
the slope correspond to lower luminosity evolution parameters and steeper 
luminosity functions (model G with cr 7 = 0.3, and model H with a 7 = 0.5). 

The best description of the observed data has been obtained by de- 
creasing at low redshift the luminosity evolution of the bright QSOs. This 
has been obtained by including a dependence on the redshift and absolute 
magnitude of the evolution parameter k such as: 
for Mb < M* : k = h + k 2 (M B - M*)e~ z / A0 
for M B > M* : k = k x 
where M* is the magnitude of the break in the two power law shape of the 
LF. This luminosity dependent luminosity evolution (LDLE) model (model 
C in Table 1) has resulted in a better fit of the data (see Fig. 3) giving a \ 2 
test probability of 0.53 in the whole B,z plane, and an acceptable 2D KS 
probability of 0.09 in the redshift interval 0.3 < z < 0.6. With this model, 
in this redshift interval and for magnitudes brighter than Mb = —25, 29 
QSOs are expected in comparison with the 32 observed. 
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Figure 5. The QSO LF. Dotted lines: best fit LDLE, model C (see Table 1). 



In this way, as shown in Fig. 1, the LDLE reproduces the higher counts 
of bright QSOs at B < 17.0, discovered by the HBQS and EQS surveys, 
much better than the PLE models. Adopting a q = 0.1 Universe, as shown 
in Table 1, the fitting the data with the PLE model (I) by Boyle (1992) is 
rejected at 0.004 confidence level. Our PLE model (L) obtains an acceptable 
global representation of the data (x 2 probability 0.36), slightly better than 
the PLE fit in a q Q = 0.5 Universe. However, in the 0.3 < z < 0.6 redshift 
range the PLE still underestimates the expected number of QSOs, although 
the discrepancy is less serious than in the q a = 0.5 case: 44 QSOs are 
observed in this domain versus 33 predicted in our L model (24 in the I 
model) . The introduction of a LDLE parameterization provides satisfactory 
representation of the observations in the low-z domain also for q = 0.1. 
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